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ABSTRACT 

The spectral distortion of the cosmic background radiation produced by the inverse 
Compton scattering on hot electrons in clusters of galaxies (thermal Sunyaev-Zel'dovich 
effect) is calculated for arbitrary optical depth and electron temperature. The distor- 
tion is found by a numerical solution of the exact Boltzmann equation for the photon 
distribution function. In the limit of small optical depth and low electron temperature 
our results confirm the previous analyses. In the opposite limits, our method is the only 
one that permits to make accurate calculations. 

Subject headings: cosmic microwave background — cosmology: theory — galaxies: 
clusters: general — methods: numerical — scattering 



1. Introduction 

The frequency spectrum of the Cosmic Microwave Background Radiation (CMBR) is known 
to have the equilibrium black body form (in natural units in which h = c = ks = 1) 

/oW = [ ex v(Pj/ T i) ~ . (!) 
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with the temperature T 7 = 2.725 ± 0.002 K (Mather et al. 1999). No deviation from this perfect 
Planck distribution is observed with the accuracy of 10 -4 . These data present a strong evidence in 
support of Big-Bang cosmology. It is established that the cosmic plasma in the early universe was 
in a thermal equilibrium state and the spectrum remained practically undisturbed to the present 
epoch. Nevertheless, small deviations from the perfect Planck spectrum are possible, and they could 
provide interesting information about physical processes in the early universe that might take place 
at relatively small red-shifts, z < 10 7 . At higher red-shifts all distortions of the thermal equilibrium 
spectrum would be efficiently smoothed down by the Compton scattering, 7 + e <-> 7 + e and by 
the inelastic photon producing reactions, double Compton, 7 + e <-> 27 + e or Bremsstrahlung, 
e + A<r+e + A + r y. The first elastic process restores kinetic equilibrium, i.e. it forces the photon 
distribution function, / T , to take the Bose-Einstein form with a possible non-zero chemical potential, 
/i, while the other two reactions push fx down to zero. 

There are several possible sources and mechanisms that could give rise to spectral distortions 
of the CMBR both in the early and the present day universe. In the early universe there could be 
electromagnetic decays of long-lived particles with life-time larger than ~ 100 sec, in particular, 
vh — ► vj_, + 7. A study of a possible distortion of the CMBR spectrum permits to obtain strong 
bounds on the probability of such decays (for a review, see Raffelt (1996)). In the present-day uni- 
verse a very interesting spectral distortion can be induced by the scattering of cosmic background 
radiation on hot electrons in galactic clusters, known as the Sunyaev-Zel'dovich (SZ) effect (Zel- 
dovich & Sunyaev 1969a,b; Sunyaev & Zeldovich 1970a,b). Observations of such a distortion, 
combined with those of thermal X-ray emission of the cluster gas, could help to extract important 
astrophysical information, in particular to measure the Hubble constant, or to study the evolution 
of clusters (for reviews see Zeldovich & Novikov (1983); Sunyaev & Zeldovich (1981); Rephaeli 
(1995a); Birkinshaw (1999a), and the updates in Birkinshaw (1999b) and Rephaeli (1999)). 

Compton scattering conserves the number of photons, which implies that the SZ effect produces 
a systematic shift of photons from the low-energy part to the Wien side of the Planckian CMBR 
spectrum. To calculate the corresponding distortion of the spectral distribution one should use the 
Boltzmann kinetic equation for the distribution function of photons, / 7 , 

f -/-, + *, (2) 

with properly taken collision integral I co u and a possible source term, S. 

In the early universe case one has to solve the system of such two coupled equations for the 
photon and electron/positron distributions. In the lowest order in the fine-structure constant, 
a ~ 1/137, only Compton and Coulomb scattering need to be taken into account. They have 
cross-sections of the order a 2 , while inelastic reactions contribute at most to the order a 3 . Because 
of technical difficulties this problem was treated only approximately especially for relativistic elec- 
trons. In the non-relativistic case the problem can be reduced to the partial differential equation 
describing diffusion in photon momentum space. For the particular case of Compton scattering 
such a reduction was done in 1957 by Kompaneets (Kompaneets 1957) and for inelastic processes 
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(double Compton scattering and Bremsstrahlung) by Lightman (Lightman 1981). Solutions to 
these equations in cosmological situations were analyzed in Bernstein & Dodelson (1990) and Hu 
& Silk (1993). 

In the case of scattering of the CMB radiation on hot electrons in galactic clusters (thermal 
SZ effect) the Kompaneets equation has the form 



is the dimensionless photon momentum, n e is the electron number density, and ctt = 87ra 2 /3m 2 = 
6.65 x 10 -25 cm 2 is the Thomson cross-section. However, the electrons in the clusters are hot, with 
temperatures that can be larger than 15 keV, which means that the energy change of Compton 
scattered photons is not small enough to be accurately described by the Kompaneets equation. 

There are basically two different approaches in the literature to extend the validity of this 
equation. In ref. Rephaeli (1995b) the relativistic form of the Maxwell velocity distribution of 
electrons is used and the frequency redistribution through Compton scattering is calculated (see also 
Wright (1979); Rephaeli k Yankovitch (1997); Sazonov & Sunyaev (1998); Molnar & Birkinshaw 
(1999)); this method is also known as the Radiative Transfer Approach. On the other hand, one 
can obtain a relativistic generalization of the Kompaneets equation, by expanding in series of the 
parameter 8 e = T e /m. Low orders of this expansion were considered in the papers Stebbins (1997) 
and Challinor &; Lasenby (1998), while Itoh et al. (1998) took into account relativistic corrections 
up to 0(#e)- These two different methods are essentially consistent. 

It has, however, been pointed out (Challinor & Lasenby 1998) that these results are not 
guaranteed to be accurate because the series approximation to the solution of the Boltzmann kinetic 
equation involves expansions in parameters which are not small. In particular it was recognized 
that the convergence is slow and possibly even asymptotic (for instance, to calculate the crossover 
frequency, where the thermal distortion vanishes, it is better to use a linear approximation in 9 e 
than the expansion up to 0(9^) or O(0^)). In view of that it is desirable to solve the original kinetic 
equation (2) directly without any specific approximation. This problem was partly addressed in 
Itoh et al. (1998, 2000a) where the collision integral was numerically calculated for the unperturbed 
photon distribution eq. (1). This integral determines the first time derivative of the distribution 
function fj(X, 0) and in the limit of a small optical depth, 



(the integral is taken along the line of sight through the cluster) gives the solution, f 7 (X, r) « 
fo + t I co u[fry(X, 0)]. In a subsequent paper (Itoh et al. 2000b), the second order corrections 




(3) 



(4) 





(5) 
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in r were calculated as a series in 9 e . The results show a good agreement with solutions of the 
generalized Kompaneets equation in the limit of small optical depth and low electron temperatures, 



In this paper we present an accurate numerical solution of the exact Boltzmann kinetic equation 
for arbitrary optical depth and electron temperature. In the limit of small r and 9 e our results 
confirm those of the previous papers, and in particular of Itoh et al. (1998, 2000b). However, the 
method presented here permits to solve the equation precisely for an arbitrary optical depth and 
any isotropic electron distribution function. We use essentially the same method as the one we 
have developed for the calculations of the spectral distortion of light (Dolgov et al. 1997, 1999a) 
or heavy (Dolgov et al. 1998, 1999b) neutrinos at Big Bang Nucleosynthesis. An essential point of 
the calculations is an analytical reduction of the exact collision integral down to two dimensions 
(Semikoz &: Tkachev 1995, 1997). In the case of Compton scattering the problem is much more 
complicated because the squared amplitude of this process is not a simple polynomial function of 
particle momenta as was the case for weak interaction in the low energy limit. 

A reduction of the collision integral down to one dimension in the direct reaction term and 
down to two dimensions in the inverse reaction term was done in Poutanen & Svensson (1996) for 
the case of Boltzmann statistics (see also the review Nagirner & Poutanen (1994) for a discussion 
and relevant references). The method of these works is somewhat different from ours because we 
specially use the procedure of the integration suitable for the quantum statistics case when the 
direct reaction term (as well as the inverse one) contains the product of the distribution functions 
not only in the initial state but also in the final state, /i/2(l + /3)(1 — fi) (see eq. 9), that makes 
the reduction down to one dimension impossible in principle if all the functions are considered as 
unknown. 

The paper is organized as follows. In the next rather technical section we present the reduction 
of the 9-dimensional collision integral down to a 2-dimensional one. In section 3 the numerical 
solution of the integro-differential Boltzmann equation is described. In section 4 our results are 
presented and discussed. Finally in section 5 we give the conclusions. 



We consider the Boltzmann equation (2) for the distribution function of photons, / 7 , taking 
into account only Compton scattering process, 



where Pi = (Ei, pj) are the particle 4-momenta. In this case the source term in the kinetic equation 
is absent, S = 0, and we have 



9 e = T e /m < 0.03 (i.e. T e < 15 keV). 



2. The Boltzmann equation 



7 (Pi) + e(P 2 ) ~ 7OP3) + e(P 4 ) • 



(6) 



#7 r 
dt ~ colU 



(7) 
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where the collision integral I co u takes the form 

Ioou = ± /n((2^) (2vrmP 1 + P 2 -P3-P4)|M| 2 F(l,2,3,4), 
with the statistical factor, F, given by 

F = / 7 (p 3 )/e(P4)[l + / 7 (Pl)][l " /e(P2)] " / 7 (Pl)/e(p 2 ) [1 + / 7 (ps)][l " AW] , 

and |M| 2 is the matrix element squared of the process in eq. (6) 

1 2 1 



(8) 



(9) 



\M\ 2 = Ae A {m A 
+ 2m 

where m is the electron mass. 



(Pi ■ P 2 



L(Pi-p 2 ) (Pi-p 4 ) 



(p 1 -p 2 )(p 1 -p 4 ) + (Pi-p 4 ) 2 
i 1 + (p • p 4 ) + (Pi • p 2 ) 



(Pi • p 2 ) (Pi • p 4 ) 



(10) 



The collision term in eq. (8) is in principle a 9-dimensional phase space integral, that can be 
analytically reduced to a sum of two-dimensional integrals. When one defines the terms 1^ for 
k = 1, . . . , 7, corresponding to the 7 terms in eq. (10), each of them reduces to 

* <-°° p 2 dp 2 rS 



r , , <** r°° p 2 dp 2 r 

h{Pl) = 7, 2 / / dP3 F J k{Pl,P2,P3) 

2vrpf Jo F 2 Jo 



(11) 



where a ~ 1/137 is the fine-structure constant and the maximum value S for the momentum j»3 
is discussed in appendix B. The expressions for the integrands Jk are listed in appendix A. Each 
of the two-dimensional integrals in eq. (11) can be found numerically, which permits to evaluate 
I co ll in eq. (8). The spectral distortion of / 7 is then found from the evolution of the Boltzmann 
equation, eq. (7) (see the next section for the computational details). An important issue, which we 
have taken into account, is the reduction of the integration region in eq. (11). When one imposes 
conservation of 3-momenta, and in particular the fact that none of the momenta can be larger than 
the sum of the others, i.e. pi < p a + p^ + p c , the integration region can be divided into several 
sub-regions where the functions Jk take different values. See the Appendices for the technical 
details. 



3. Numerical solution of the Boltzmann equation 

We have found the spectral distortion of / 7 produced by the SZ effect from the numerical 
solution of the Boltzmann equation (7). 

We measure all momenta in units of photon temperature. Accordingly, we use the dimensionless 
photon frequency X from eq. (4) and measure time in units of the optical depth, r, given by eq. (5). 
In order to calculate the integral over electron momenta we introduce the electron temperature as 

Xe = = RT^ , (12) 
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where 

Rt = ^ , (13) 

is the ratio of electron to photon temperatures. We assume that electrons are always in equilibrium 
with the distribution 

f e = (e^e-^/^ + l) -1 , (14) 

where E e = ^/m 2 + i s the electron energy and fi e is their chemical potential. In the present 
paper we make the simplifying assumption of Boltzmann statistics 4 for electrons (f e <C 1) and 
express fi e through the number density of electrons. 

The integration over electron momenta in the two-dimensional integral in eq. (11) should 
formally be taken from zero to infinity, but we integrate only to p™ ax . We found that for small 
T e < 10 keV it is enough to take p™ ax = m, and for larger electron temperature T e ~ 50 keV 
integration till p™ ax = 3m is sufficient. 

For the dimensionless photon frequency, X, we introduce a grid of iV gr id points, logarithmically 
spaced in the interval X\ ow < X < X^. We use X\ ow = 10 -4 , which is small enough to maintain a 
sufficiently high precision of the calculations. X^ actually depends both on the electron temperature 
T e and the final optical depth Tfi n . For small electron temperatures T e < 10 keV it is enough to 
take Xhi = 100. For higher temperatures and large optical depth r > 1 we need to use much 
higher values of X^, because photons multiple-scatter into large momentum modes. For example, 
for calculation with T e = 20 keV and Tfi n = 20 (see the next section) we use X^\ = 10000. As for 
the number of points in the grid, iVgrid, we always checked that it is large enough to maintain the 
desired precision. Typically 1600 points grid was good enough for our purposes. 

As initial condition for eq. (7) we used equilibrium photon distribution function eq. (1), which 
in terms of X has the simple form 

MX) = . (15) 

For calculation of the two-dimensional integrals in eq. (11) we use the Gaussian quadratures 
method, see Press et al. (1992). For the evolution in "time", r, we use the simple Euler method, 
which is precise enough for the present calculations. Finally, for calculation of the photon distri- 
bution function in points between the grid points we used interpolation not in / 7 but in / 7 — /o to 
perform 2D-integration. This trick strongly reduces numerical errors, coming from interpolation in 
the regions where the distribution of photons is close to the equilibrium one (i.e. for r < 1). 



4 Note that our method is valid for any isotropic / e , not necessarily in equilibrium, but here we restrict ourselves 
to the case of a Boltzmann equilibrium distribution in order to investigate the thermal SZ effect. 
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Additional difficulties in the numerical calculations come from the region of small photon 
momenta, X < 0.1. In this region some parts of the integrands in eq. (11) behave like 1/X 5 , while 
the complete collision integral is proportional to 1/X for small X (due to the singularity l/X in 
the distribution function eq. (15)). This leads to large numerical errors for small momenta. On 
the other hand, asymptotically for small X each part of the collision integral eq. (11) has either 
a constant value or 1/X-behavior. Using this fact, we calculate the integrals numerically for grid 
points with momenta X > X smaj \\ and extrapolate into the region of smaller X. In the calculations 
we use 0.01 < X sma \\ < 0.1. 

Different parts of the integration region contribute differently to the collision integral in eq. (11) 
(see appendix B for the definitions of the regions). The main contribution comes from the large 
regions C2 and Z?2- The regions C\ and D\ can be taken into account, but they contribute on the 
level of 1% in the region of small X < 0.1 and could be neglected for larger X. The main numerical 
error comes from a small part of the integration regions C and D around the line Q 2 = 0. These 
regions give numerical errors up to 10% of the value of the collision integral for small photons 
momenta. The nature of this small-momentum error was discussed above. The regions A and B 
give corrections to the collision integrals in the sixth digit and could be neglected. 

The ratio of the electron and photon temperatures, Rt (13), is very large in hot clusters. 
The CMB temperature is T 7 ~ 10 -4 eV, while the electron temperature is T e ~ 10 keV. Since 
both photon and electron momenta are involved in the same expressions in the collision integral 
we may lose precision by many orders of magnitude if we keep realistic values of Rt ~ 10 8 . The 
energy conservation implies p% = p\ + E2 — E4, and since E2 ~ E4 ~ m = 511 keV, the difference 
should cancel up to 10 digits in order to give numbers of the order p\ ~ p% ~ 10~ 4 eV, and even 
more precisely for low photon momenta, X < 1. This numerical problem is solved in the following 
way. We calculate the collision integral for much higher photon temperatures T 7 ~ 1 keV, and 
then diminish the temperature until the collision integral stops varying with a further temperature 
decrease. We assume that at this value of T™ 111 the collision integral reaches its asymptotical value 
and therefore reproduces the result for T 7 ~ 10~ 4 eV. We found that already at T 7 = 50 eV 
the collision integral reaches asymptotical values. In most calculations we used T™ m = 10 eV. 
Let us note, that already for T 7 = 1 eV the numerical errors discussed in this paragraph become 
important. 

In conclusion, we control all numerical errors in the solution of the kinetic equation (7), and 
we also check the result testing the conservation of the photon number (see the end of section 4) . 

4. Results 



The distortion of the photon spectrum is usually presented in the form 

AI(X,T) = X 3 (f^X,T)-MX)) , 



(16) 
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where fy(X, r) is the photon distribution function calculated from the kinetic equation (2) and 
fo(X) is the initial equilibrium photon distribution (1). 

Instead of solving the kinetic equation directly, one could for a small optical depth, r <C 1, 
make the Taylor expansion in r of the collision integral 

I coll (X, t) « 7 coH (X, 0) + r^(X, 0) + 0(r 2 ) . (17) 

Then the distribution function fj(X, r) can be trivially found from eq. (7) 

/ 7 (X,r) « /o(l) + rl mi ,(l,0) + ^^(X,0) + 0(r 3 ) , (18) 
and the spectral distortion takes the form 

AI(X, t) « A/i(X, r) + A/ 2 (X, r) + 0(t 3 ) , (19) 

where 

AI 1 (X,t)=tX 3 I coU (X,0) , A/ 2 (I,T) = i 3 ^(I,0) (20) 



The term AIi in eq. (20) corresponds to the single photon scattering. It requires only calcula- 
tion of the collision integral with equilibrium distribution functions fo(X). Our calculation of the 
collision integral in this case agrees perfectly with the results obtained in Itoh et al. (1998, 2000a), 
as can be seen from fig. 1. In this figure we compare our numerical results for T e = 10 keV and 
T e = 20 keV with analytical fits from Itoh et al. (2000a). 

The term A/2 in eq. (20) is proportional to r 2 , and it corresponds to the double photon 
scattering contributing to the spectral distortion. In order to calculate this term one can do the 
following: calculate I co u(X, 0), substitute the first two terms of the photon distribution function 
eq. (18) back into the collision integral and calculate all terms of order r. These terms correspond to 
the second term in eq. (17), which is responsible for the double photon scattering. In the paper Itoh 
et al. (2000b) this second order contribution to the spectral distortion was found approximately, up 
to the terms of the order Q 2 e . By numerical integration of eq. (7) we have found the two-scattering 
contribution exactly from I co u (X, r) . 

In fig. 2 the double-scattering contribution into the spectral distortion as a function of dimen- 
sionless photon frequency X is presented for T e = 10 keV. The solid curve in fig. 2 is our exact 
result (up to small numerical errors). One could compare it with the series expansion in powers 
of 9 e = T e /m from Itoh et al. (2000b), which are shown by the dashed, short dashed and dotted 
lines (the definitions of functions Z{ are given in Itoh et al. (2000b)). As we can see from fig. 2, 
this expansion is in a good agreement with our result for relatively small X < 4 and disagrees for a 
larger X. It demonstrates again (as the first-order terms in r) a bad convergence of the expansion 
in terms of 8 e . 
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In figure 3 we present the results for the spectral distortion of the CMB photons caused by 
the SZ effect for T e = 20 keV and the optical depth equal to one, r = 1. The dotted curve 
corresponds to the numerical solution of eq. (7). The solid line represents the single scattering 
contribution, namely the first term in eq. (19). The long dashed line corresponds to the double 
scattering contribution, namely the second term in eq. (19). The short dashed line is the sum of 
these two contributions. One can see from fig. 3 that even at optical depth r = 1, the spectral 
distortion is accurately described by the sum of single and double photon scatterings. 

It should be noted that our results have been obtained with an initial isotropic photon distri- 
bution. This is an approximation of the real radiation distribution inside the clusters of galaxies, 
which are of finite extension. At the level of single scattering there is no problem since the radiation 
is isotropic prior to convolving with the Compton kernel. However the multiple scattering contri- 
butions to the spectral distortion could be affected by the geometrical properties of the cluster. 
For instance, Molnar & Birkinshaw (1999) used a Monte-Carlo method to take into account finite 
optical depth and bulk motion in a spherically symmetric uniform plasma. Therefore our results 
are exact in the case of an infinite medium. 

For r <C 1, the single scattering strongly dominates the spectral distortion, and in order to 
estimate the contribution of the second-order term, one should rescale the amplitude of the long 
dashed line in fig. 3 by the factor r and compare it with the solid line at any photon frequency X. 
In a realistic case one has r < 0.01, and the double scattering can be safely neglected everywhere, 
except for a small region in photon frequency X where the single scattering contribution vanishes. 

The magnitude of the double photon scattering is demonstrated in figure 4 for the optical 
depth t = 0.01. The crossover frequency Xq (the photon frequency at which the distortion caused 
by the thermal SZ effect vanishes) is plotted as a function of the electron temperature (solid line). 
The regions around this line show the double scattering contribution to the spectral distortion. In 
the region between the dotted lines it is larger than 50%, between the short dashed lines it is larger 
than 10%, and between the long dashed lines it is larger than 1%. 

When one calculates the distortion using the Kompaneets equation (3), the crossover frequency 
does not depend on the electron temperature, and it takes the constant value Xq ~ 3.83. The 
exact calculations show that the crossover frequency depends both on the electron temperature, 
6 e = T e /m and the optical depth r. For small optical depths, < r < 0.05, a linear fit in r is very 
accurate 

X = a(T e ) + t (3{T e ) , (21) 

and by doing such fits in the range, < r < 0.05, for many different electron temperatures one can 
find a(T e ) and /3(T e ) as functions of temperature 

a(T e ) « 3.830(1 + 1.162 9 e - 0.8144 6 2 e ) and (3{T e ) « 3.021 6 e - 8.672 6 2 e , (22) 

which fits better than 4 x 10~ 4 for < T e < 50 keV and < r < 0.05. The functional shape 
was chosen by simply looking at the corresponding graphs. By considering the graphs for a(T e ) 
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and P(T e ), we have found no reason to include higher order terms. The fit parameters change only 
slightly for large optical depths, 0.05 < r < 1. One can naturally perform similar fits for X max and 
-Xmin) which are the dimensionless frequencies where the spectral distortion is maximal or minimal. 
We find 

X min = 2.265 (l- 0.0927# e + 2.38 0*) + r (-0.00674 + 0.466 e ) , (23) 
X max = 6.511 (l + 2.41 9 e - 4.96 0j) +r (0.0161+ 8.16 e - 35.9 0j) . (24) 

We can compare these results with similar fits from the previous papers. For the interesting range 
of temperatures, < T e < 50 keV, and for negligible optical depth, we are in good agreement with 
the previous results (Birkinshaw 1999a; Challinor & Lasenby 1998; Itoh et al. 1998). In Molnar & 
Birkinshaw (1999) the fit in the optical depth is also made and for small r we are in good agreement. 
In this paper some of the effects of the cluster geometry were considered, which may explained why 
our results for larger r are somewhat different. 

Finally in figure 5 the results for a large optical depth, r, are presented. The Taylor expansion 
in terms of r, eq. (19), is useless for r > 1, and the complete numerical solution of eq. (7) should 
be done. From figure 5 one can see that for r > 1 the spectral distortion is dominated by multi- 
scattering contributions. Indeed, one and two-scattering contributions dominate the curves with 
r = 0.01 and r = 1, and all distortion in these cases are in the region X < 20. Spectral distortion 
in higher modes come from multi-scattering contribution, demanding much larger optical depth. 

Asymptotically for r — > oo the photon distribution should reach the form 

(25) 

with the chemical potential p determined by the photon number conservation, p = T e log(C(3)T^/T^) ~ 
— 1 MeV (when T 7 <C T e ). However, even for r = 20 the distribution is very far from the equilibrium 
one in eq. (25). A much larger r is necessary to equalize the temperatures of photons, T 7 = 2 x 10~ 4 
eV, and electrons, T e ~ 20 keV. We have estimated that one requires r ~ ^m/T e ln(p ^ n / pi n ) , 
where pi n and pfi n are the initial and final energy densities of photons. We have also calculated 
the same quantity numerically and found r ~ 250 for T e = 20 keV, in reasonable agreement with 
this simple analytical estimate. As expected we observe that the crossover frequency increases with 
rising r, and asymptotically it should tend to infinity. 

To check the precision of the calculations we test if the number density of the photons is 
conserved, as it should be in elastic Compton scattering. For T e = 20 keV and 1600 point grid in 
photon frequency the relative non-conservation of the photon number, AN^/N-y, is approximately 
10~ 9 for r = 0.01 and is smaller than 10 -5 for any r < 20. 
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5. 



Conclusion 



We have numerically calculated the spectral distortion of the cosmic background radiation 
produced by Compton scattering on hot electrons in galaxy clusters without any simplifying as- 
sumptions. We have analytically reduced the collision term in the Boltzmann equation to two 
dimensional integrals, and accurately solved the kinetic equation. Our analytical expression for the 
collision integral is exact, and can be used for any temperatures, T 7 , T e , and optical depth, r. The 
results are in good agreement with the previous works for small temperature and optical depth, 
but considerably different for large r. Our method can also be applied in the case of an isotropic 
electron distribution function, not necessarily in equilibrium. 
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by INFN. In Munich, this work was partly supported by the Deutsche Forschungsgemeinschaft 
under grant No. SFB 375. The work of D.S. is supported in part by INTAS grant 1A-1065. S.P. is 
supported by the European Commission under the TMR network grant ERBFMRX-CT96-0090 
and a Marie Curie fellowship under contract HPMFCT-2000-00445. 



In this appendix we list the expressions for the 7 terms in the collision integral corresponding 
to the 7 terms in the matrix element eq. (10) 



The functions J±, . . . , J7 depend upon the following combinations of the particle 3-momenta 



A. Reduction of the collision integral 




(Al) 



V 2 = (pi + E 2 f - m 2 
a = max(|pi - p 2 \, \p3 ~ Pi\) 
c = max(|pi -p 4 |, \p 2 — P3I) 



Q 2 = (E 2 - p 3 ) 2 - m 2 , 
b = min(pi + p 2 ,P3 + Pa) , 
d = min(pi +p±,p 2 + Pz) ■ 



(A2) 



In each case the value of p^ is obtained from the condition of energy conservation 




(A3) 



The integration over d?p^ and the angles of P2 and P3 with <5-function imposing energy-momentum 
conservation is straightforward but rather lengthy, and we present only the final results 



T 4 

J\ = m 




(A4) 
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(A \ (h \ 

J 2 = (d-c) = (b - a) , 



J 4 



PlPZ 

2m 2 



V 

b — a 



log 



PlP3 

(V + b){V-a) 



1 + 



(V-b)(V + a 

(p? -pI)(pI-pD 



V 2 ab 



m 2 p 1 p 3 
+ —ps— } °Z 



(V + b)(V-a) 



(V-b)(V + a) 



The explicit form of the other three Jj depends on the sign of Q 2 . If Q 2 > then 



T 4 

Jo = m 



0^ 1Og 



(Q + d)(Q-c) 



2m 2 
d — c 



log 



(Q-d)(Q + c) 
(Q + d)(Q-c, 



+ 



2(d - c) (1 + cd/Q 2 
(Q 2 -d 2 )(Q 2 -c 2 ; 



(Q-d)(Q + c) 



1 + 



+ 



2 \ Q 2 c(i / Q 3 

while for the case Q 2 < 0, defining W 2 = -Q 2 we find 



log 



(Q + d)(Q-c) 



{Q-d)(Q + c) 



h 
J 7 



2rrf 
4m 2 



_L arctan + (^-c)(i + ^/Q 2 ) 



d — c 



arctan 



(d - c)W 



1 + 



cd- Q 2 
Q 2 cd 



2m z pipz (id - c)W 

■ arctan 



Q 2 W 



cd - Q 2 



(A5) 
(A6) 
(A7) 



(A8) 

(A9) 
(A10) 



(All) 

(A12) 
(A13) 



B. Integration region 

In figs. 6 and 7 we present the integration region for the two-dimensional integral (Al), where 
we took the particular case p\ = l,m = 50 p\ and measure momenta pi in units of p±. The line 
j>3 = S = p\ + E% — m is the formal upper limit of integration over p^, so < p% < S. The 
integration region for p2 is < P2 < oo, but in fig. 7 we present it only for p2 < m. Throughout 
the paper we consider the case of non-relativistic electrons, p e < m but since we never expand in 
the electron mass, our analysis is valid for the general (relativistic) case. 

We can cut away most of the integration region for p% from the condition pi < p a + Pb + Pc, 
which follows from the momentum conservation, i.e. one momentum cannot be larger than the sum 
of all the others. For pi = p\ or p?. this condition gives no reduction, but in case of pi = p2 or p4 
we cut regions E and F from the integration region 5 (see figs. 6 and 7). Thus only the regions A, 
B, C and D survive. To quantify this reduction let us note, that in the case of CMB distortion 



5 Note that region E only exists if pi < m/2. 
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in hot galaxy clusters one has p e ~ ^2mT e ~ 100 keV, while p 7 ~ 10 -4 eV. In this case formally 
< Pz < S ~ 100 eV, while non-zero contributions to the integrals in eq. (Al) come only from 
P3 ~ 10~ 4 eV. 

Figure 6 shows the integration region for relatively small momenta. One can see that besides 
the main integration regions C and D there are two small regions A and B which appear due to 
the fact that for a small momentum p 2 , the condition P2 = Pi + P3 + Pi is not operative, and the 
upper limit for p 3 becomes p% = S. 

The regions A,B,C and D are defined in such a way, that in each of them our parameters a, b, 
c and d take definite values. This allows one to simplify the collision integral, I co u, hi each region. 
As we have seen in the previous section, there are generally two kinds of contribution to the collision 
integral, like Ji (which we will call P-kind) and like J3 (which we will call Q-kind). Fortunately 
the integration regions are identical for all integrals, and only the limits a, b, c and d will differ for 
the two kinds. For Q-kind integrals we must distinguish between positive and negative Q 2 in the 
regions C and D (because Q 2 is always negative in regions A and B). Therefore we introduce the 
regions C\ and D\ for Q 2 < 0, and C2 and D2 regions for Q 2 > 0. We have also separated the 
contribution from the Q 2 = case, but this one-dimensional integral gives negligible contribution 
to the two-dimensional collision integral. Below we define each of the regions in detail. 

Region A. 

< p 2 < pi ; pi < P3 < pi + E 2 - m , 
and 



If p\ < m/2 then, 

Pi(E2 + P2) „ 

Pl<P2< Pchange 5 7, ~~ET < P3 < Pi + E 2 ~ m , 

*Pl — P2 + &2 

where p^ange = 2pi(m — pi)/(m — 2pi) is the point in which p2 = Pi + P3 and p^ = 0. 
If pi > m/2 then, 

Pi(E2+p 2 ) ^ 

pi < p 2 < oo ; — — < p 3 < pi + E 2 - m , 

2pi - P2 + E 2 

and region E does not exist. 



In region A we have 

a = P3~P4 ; b = p 3 +p i ; c = pi - p 4 ; d = pi+p A . 



Region B. 
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In this region we have 

a = Pl -p 2 ; b = pi + p 2 ; c = p 3 - p 2 ; d = p 2 + p 3 • 



Region Ci. 

In this region Q 2 < 0. It is defined as 

. ■ , n Pl(E 2 -p2) Pi(E 2 +p 2 ) 
< p 2 < min(pi, pcj 5 7; ; —ft < P3 < 



2pi + P2 + ^2 2pi - p 2 + E 2 

and 

If pi < to/2 then, 

Pi{E 2 -p 2 ) 

Pl<P2< PCi ; o T r~ET <P3<Pl, 

2pi + P2 + -&2 

If p! > m/2 then, 

Pi(£ 2 + P2) 



Pd <P2 <Pi ; E 2 -m<p 3 < 



where 



2pi - P2 + E 2 ' 
and 

max(pi,pci) < P2 < PC 2 ; E 2 -m<p 3 <p 1 , 



(m-2 Pl ) ^/pi (2m+pi ) -2pi (m-pi ) / m 

m— 4m "1 ' 



m— 4pi ~" 7^ 4 

P< : { (Bl) 

n m pi = f . 

is the point, in which the line Q 2 = crosses the line P4 = pi + p 2 + P3 (see figures 6 and 7). Let 
us just note the limiting cases for pc x 

Pi < to p Cl -> y^piTO , 

TO TO 

Pi = y PCi = Pi = y . 

3 

Pi -> 00 p Cl = -to . (B2) 

The point in which the line Q 2 = crosses the line P3 = pi is 

pc 2 = Vpi(2toTp0 • (B3) 



-15- 



In region C\ we have 

a = pi-p 3 ; b = pi + p 2 ; c = p 4 - p\ ; d = p 2 +p 3 . 



Region C 2 . 

In this region Q 2 > 0. It is defined as 



Pi(E 2 — P2) _ 
Pd <P2 < PC 2 5 — — 1 1 1 <Ps < E 2 -m, 



and 



Pc 2 < < 00 



2pi + P2 + -^2 



2pi + j> 2 + £2 
The limits a, b, c and d are the same as in region Ci. 



< P3 < Pi 



Region Di. 

In this region Q 2 < 0. It is defined as 



Pi <P2 < Pc 2 5 Pi < Ps < 



Pi(E 2 +P2) 



and 



PC 2 < Pi < PD ; E 2 -m < p 3 < 



Pi(-E 2 + P2) 



2pi - P2 + -^2 ' 

where po is the point where the Q 2 = line crosses the line p 2 = P1+P3+P4. This point is defined 

as 



PD = < 



(m-2pi ) -y/pi (2m+pi ) +2pi (m— pi ) 
m— 4pi 

3m/(4e) 



00 



Pi < 



Pi = x - me , e < 1 



(B4) 



In region Di we have 

a = p 2 - pi ; b = p 3 + p 4 ; c = p 2 - P3 ; d = p x + p 4 . 

Region D 2 . 

In this region Q 2 > 0. It is defined as 

Pc 2 < P2 < PD ; pi <P3 < E 2 -m, 
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and if po is finite 

Pl(^2+P2) 

£>D < £>2 < °o ; pi < _p 3 < 



2pi - P2 + #2 
The limits a, b, c and d are the same as in D\. 
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Fig. 1. — First-order contribution to the spectral distortion of photons for T e = 10 keV and T e = 20 
keV as a function of the dimensionless photon frequency X. Our numerical results are presented 
by the solid and dashed lines, whereas the data points correspond to the analytical fit from Itoh et 
al. (2000a). 

Fig. 2. — Second-order contribution to the spectral distortion of photons for T e = 10 keV as a 
function of the dimensionless photon frequency X. Our numerical result are presented by the solid 
line, whereas the dashed, short dashed and dotted lines correspond to the series expansion from 
Itoh et al. (2000b). 

Fig. 3. — Spectral distortion of photons for optical depth r = 1 and T e = 20 keV as a function of 
the dimensionless photon frequency X. 

Fig. 4. — Crossover frequency Xq as a function of the electron temperature T e for r = 0.01 and the 
regions, where double photon scattering gives a contribution to the spectral distortion larger than 
1% (long dashed lines), 10% (short dashed lines), and 50% (dotted lines) respectively. 

Fig. 5. — Spectral distortion of photons for T e = 20 keV as a function of the dimensionless photon 
frequency X for several optical depths up to r = 20. 

Fig. 6. — Region of integration of two-dimensional integral for small momenta. 
Fig. 7. — Region of integration of two-dimensional integral. 
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